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We explore the sensitivity of weak gravitational lensing to second-order corrections to the space- 
time metric within a cosmological adaptation of the parameterized post-Newtonian framework. 
Whereas one might expect nonlinearities of the gravitational field to introduce non-Gaussianity into 
the statistics of the lensing convergence field, we show that such corrections are actually always 
small within a broad class of scalar-tensor theories of gravity. We show this by first computing the 
weak lensing convergence within our parameterized framework to second order in the gravitational 
potential, and then computing the relevant post-Newtonian parameters for scalar-tensor gravity 
theories. In doing so we show that this potential systematic factor is generically negligible, thus 
clearing the way for weak lensing to provide a direct tracer of mass on cosmological scales for a wide 
class of gravity theories despite uncertainties in the precise nature of the departures from general 
relativity. 
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I. INTRODUCTION 



Weak gravitational lensing, whereby galaxy images are altered due to the gravitational influence of the mass along 
the line of sight, is a powerful probe of the "dark sector" of cosmology, with promising results in recent years, e.g. [H- 
[S]. (See Refs. [BHH] for recent reviews.) In general relativity (GR), the weak lensing distortion field provides a direct 
tracer of the underlying matter distribution. However, modifications to gravity theory can conceivably alter the way 
that mass curves spacetime and thus the way that null geodesies correspond to the matter distribution. The linear 
order effect, often referred to as the "gravitational slip," has already been explored in the literature [9], with weak 
lensing and other datasets. But what happens at higher order? Higher-order lensing statistics are useful for breaking 
degeneracies and probing primordial non-Gaussianity (e.g. [51 IIP)), but could modifications to GR complicate such 
efforts by further altering the way that lensing traces mass? On the other hand, could we conceivably use higher-order 
weak lensing statistics to probe nonlinear gravitational dynamics? 

We aim to answer these questions by utilizing a cosmological adaptation of the parameterized post-Newtonian 
(PPN) framework |TTJ[T2]. For example, consider the following special case for the static weak-field metric: 

1 - |C/ + f C/^) c-de + (l + + I^C/^) dx-^ , (1.1) 

where U is the usual Newtonian gravitational potential and the powers of 1/c^ label the order of the perturbation 
from Minkowski spacetime. The parameters take on the values ^ — P — e = \ \n the case of GR. At second 
order, /3 and e represent the first nonlinear effects of gravity^. In cosmology we would like to determine if new laws 
of gravitation are responsible for cosmic acceleration, possibly replacing the need for dark energy 1131 114j. If new 
gravitational phenomena are at play, then we should expect a departure from GR in the cosmological analogues of 
7,/3, and e, which may be manifest in the patterns of gravitational lensing of light from galaxies, clusters, and the 
cosmic microwave background (CMB). 

The nonlinearity of the gravitational field introduced by the /3 and e terms translates into a non-Gaussianity of the 
statistics of cosmological fluctuations. Indeed, standard single-field inflation models predict nearly Gaussian initial 
cosmological fluctuations [TS] . As such, probes of non-Gaussianity provide a crucial test of the inflationary paradigm, 
and thus much effort has been put into constraining non-Gaussianity in both the early- and late-time universe. Local 




^ Because non-relativistic objects move with \dx\ <^ cdt, the e term enters at a higher order than does the 13 term. For this reason the e 
term is often neglected in investigations of alternative theories of gravity. However, the e term is necessary for a self-consistent analysis 
of the equations of motion of light, i.e. gravitational lensing. 
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type non-Gaussianity is typically parameterized by the parameter /nl, where on sub- horizon scales the Newtonian 
potential (equal to U above and below) can be re- written in terms of a Gaussian field (j)o: 



4>{x) = Mx) + /nl [M^f - {4>l{x))] 

If we take the GR metric and replace U ^ U + A/nlI^^, then we find 



1-^(7+^(1-A/nl)(72 



1 



4c/+4('- + A/nl)C/^ 



dx' 



(1.2) 



(1.3) 



to second order in U. Comparing this with the line element (1.1 1, we may infer that additional non-Gaussianity 
beyond GR is introduced by the /3 and e terms. Such additional non-Gaussianity could conceivably be confused with 
that of primordial origin. In this paper we determine how gravitational lensing can probe this effect; by computing /3 
and £ for a wide class of scalar-tensor gravity theories, we further show that such corrections to /nl are generically 
negligible for the majority of the modified gravity literature. In doing so, we further show how the second-order 
lensing convergence depends on the "gravitational slip" 7. 

The organization and results of this paper are as follows. In Section II, we compute the weak lensing convergence 
to second order in the gravitational potential U, or to 0{l/c*), within a cosmological extension of the PPN formalism 
so that we may find how nonlinear lensing and non-Gaussianity depend on modifications of GR. Our first main result 
is our final answer for the weak lensing convergence, which is is given in Eqs. ( 2.26||2.29 ). Note that throughout this 
paper 7 will refer to the PPN parameter and not to the cosmic shear field. Then in Section III we study a wide class 
of scalar-tensor theories of gravity, wherein we provide a derivation of the relevant post-GR parameters and show 
that the combination of parameters important for lensing negligibly differs from the GR expectation. We conclude in 
Section IV with a discussion of our results. 



II. NONLINEAR GRAVITATIONAL LENSING 



Generalizing the PPN metric (|1.1|) to allow for non-static gravitational potentials, we write 
2 



ds' 
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2/3 



$2 c^dr 



r Vi cdt dx^ 
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27 



U 



3e 
2? 



$2 dx' 



(2.1) 



to order 0{l/c^) [HI [T^], where the potentials are defined such that GR is recovered when j = f3 = e = f3' = e' = 1, 
the necessary PPN parameters are absorbed into the definition of the gravitomagnetic potential Vi, and $2 is the 
post-Newtonian scalar potential as in |12j . To start, we assume that the background spacetime is Minkowski, but we 
will generalize to the case of an expanding Universe at the end of the calculation. We also now set G = c = 1, but 
leave the powers of 1/c in place for post-Newtonian bookkeeping purposes. 
The propagation of photons is determined by the null geodesic equation 



dfc" 
~dX 



— ipyft- ft. 



(2.2) 



where A:" (A) = dx°'/dX is the photon 4-momentum, a;" (A) is the light ray trajectory, A is an affine parameter, and 
the connection is for the full metric. Let us assume that the unperturbed ray moves along the positive x-axis, and we 
will place the source at a; = and the observer at some location x, corresponding to a comoving distance w{x). We 
enforce the null geodesic condition, ^^j/fc^fc" — 0, and normalize = (fc*, k^, k^, k^) = (1, 1, 0, 0) for the unperturbed 
ray. The geometry of the generic lensing scenario is sketched in Figure [T] The angle a is the deflection angle felt by 
the ray and a is the angular change in sky position as seen by the observer. 

We can change the independent variable from A to 2: using the following relation: 



Thus, if we define 



d dx d d 
dX dX dx dx 



W''{x) = -- 



(2.3) 



(2.4) 



as a function of the independent variable x, 



then we can integrate Eq. (|2.2|) along the ray trajectory, with the help of Eqs. ( 2.3 1 and ( 2.4 ), to find the 4-momentum 

k°'{x) = / dx'W^lxix')] . (2.5) 
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Observer 




Source 

FIG. 1: A sketch of our setup. The observer is a comoving distance w away from the source, 9 is the angular coordinate on 
the sky, a is the deflection angle felt by the ray, and a is the angular change in sky position as seen by the observer. 



Here we are integrating along the perturbed ray trajectory x(x), 

which we are also parameterizing as a function of the x-coordinate along the ray. Within the perturbative post- 
Newtonian framework, these equations can be solved iteratively to get the full ray trajectory x{x) to any order in 1/c. 
Once we know the ray position at the observer, x^{x), then the angle a* (see Figure [l]) is such that x^{x) = xa^{x), 
where now i = y, z ii we use the flat sky approximation. Thus 

"^(^) = Z r rllm r dx"W^ [x{x")] (2.7) 



X Jo [x {x')] Jo 

is the deflection angle seen by the observer. 



A. The deflection angle 

Generalizing to the cosmological scenario, whereupon the metric is replaced by 5^1/ — )■ a?'{t)g^^ so that all times 
are conformal and coordinate positions are comoving, we now denote w as the comoving distance to the image source 
and 9 as the location on the observer's image plane. Then we find 

a.(^,w) = — (1 + 7) / dw' ) d^U\w'B, w'\ + 2PN terms + . . . (2.8) 

c Jo V ^ / 

to first post-Newtonian (IPN) order. This is the standard result [16] . Note that we have pulled 7 out of the integral, 
thus ignoring its possible spatial dependence. Generally 7 will take on different values on sub-galactic scales than on 
larger scales so as to maintain the GR limit in the well-tested very small-scale regime. However, such a small-scale 
transition can be safely ignored for the purposes of the weak lensing observations discussed here, which are limited in 
resolution to scales of ~ 1 arcminute, or 0.5 Mpc at z = 1. The expression above also uses the Born approximation, 
whereby the photon moves along the unperturbed geodesic, XBorn = wQ. In fact, the linearly-perturbed geodesic is 
a^Linear = 2?Born + ^x where 5xi = —wai- To implement this correction, we Taylor expand the potential: 




+2PN terms + ... . (2.9) 
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Now letting the unit 3-vector ft denote the direction of propagation of the unperturbed ray, the displacement vector 



ai{9,w) — / dw' 
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l+7^,,r..y.- ..n (1+7)' 



-d,U[w'0,w' 



d,jU[w'0,w'] / dw"{w' -w")djU[w"9,w^ 



(^n X |v X V[w'0,w']^^ 



6 - 4/3 + 3e - 672 



2c4 



U[w'e,w']d,U[w'9,w' 



e' -p' 



d,<^2[w'e,w'] 



dw' 

w 



^-^u[w'e,w'] / dw"d,u[w"e,w" 



+ 27(1 + 7)n- V[/[w'^,w'] / dw"d,U[w"e,W 



(2.10) 



to 2PN order, where all of these potentials are evaluated along the unperturbed ray trajectory. The first line contains 
the linear order answer plus the Born correction and the second line contains the gravitomagnetic term. The first line 
is studied in Ref. [T7] and the gravitomagnetic term is studied in Ref. [18] . 

The integrands are evaluated along the line of sight, and so some simplifications may be made. Note that 



aw' 



(2.11) 



where the overdot denotes a derivative with respect to t[w') — w — w' . Then the second order term on the fourth 
line of Eq. (2.10) can be rearranged thusly: 



dw' 
dw' 
dw' 

V 
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w ~ w' 
w 

w — w' 
w 
w — w 
w 



VU[w'e,w'] / dw"d^U[w"e,w"] 



u[w'e\w'] + J^u[w'e]w'] ] I dw"d,u\w"e 

OW' j ./n 

1 



U[w'9,w']J dw"diU[w"e,w"]+ J dw' i^-]U[w'9,w'] I dw"diU[w"9,w"] 
U[w'9,w']d.,U[w"9,w"] . 



(2.12) 



Using this substitution, we find the full deflection angle to second post-Newtonian (2PN) order to be 



ai{9,w) 
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-^{l+"f)d,U[w'9,w'\ 



{l + -iYd^jU[w'9,w'] / dw" {w' -w")djU[w"0, 
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(6 - 4/3 + 47 + 3e - 2-i^)U[w'9, w']d,U[w'9, w'] 
^-f{l + j)U[w'9,w'] / dw" diU[w"9,w"] 
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+ + j^)U[w'9,w'] / dw"d,U[w"9, 



|n X (^V X V[w'9,w']^Y 



— {s'-P')d,<f2[w'9,w'\ 
\v,[w'9,w'] . 



(2.13) 
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B. An aside: Point mass case 



We can specialize the above results to the well-studied case of a static point mass. A point mass located a distance 
b away from the unperturbed ray trajectory produces a gravitational potential 

U(x, y, z) = , (2.14) 

where we have arbitrarily placed the mass at {x,y,z) = (a;;, 6,0), meaning that there will be deflection in the y- 
direction only. For the point mass, we also have 4>2 = because it is static. We further assume that the source 
and observer are very far away from the lensing mass, meaning that we can safely take the lens position to be such 
that xi ^ and the observer position to be such that a: — >■ oo. Plugging this all into the procedure outlined above 
simplifies the equations greatly, and we find the deflection angle to be 

..a.„S.(2.2.-..?.)5^ (2..) 

to 2PN order, where we have now calculated a instead of a so as to compare directly with the literature. The IPN 
piece of this answer is the widely accepted result, and the 2PN piece agrees with Ref. (TS]. 



C. Order of magnitude scalings 



A close inspection of the deflection angle formula above, Eq. (2.13) 
various terms as in Ref. [20 



permits us to estimate the scalings of the 
and determine the relative contributions. Note that w ^ th, where rn ^ 3000ft-^^ Mpc 



rn, wncre th 

is the Hubble radius, and spatial derivatives in the plane of the sky introduce factors of 1/A, where A is the mode 
wavelength. For the perturbations which contribute the most to deflections, A ~ Amax ^ 30/i~^ Mpc ~ r/f/lOO. This 
also justifies our discard of time-derivative terms, under the implicit assumption that the PPN-like parameters evolve 
on a time scale comparable to the Hubble time. Also note that the rms amplitude of the gravitational potential on 
these scales is C/rms ~ 10~^. With this knowledge in hand, Eq. (2.13) contains the following terms: 



• Line 1: j^' dw' + i)d,U[w'lw'] 

This is the leading order solution, which scales like [rn /X)U and is thus of order lO^'^ for the dominant mode. 

• Line 2: dw' (^^) ^(1 + "ffd,jU[w'e, w'] J^' dw" {w' - w")djU[w"9, w"] 

This is the Born correction, which scales like {r h / \YU'^ ■ That means that it is of order 10^^ for the dominant 
mode, making it roughly 10% of the leading order term for that mode. However, this term is significantly 
suppressed on larger scales. 



• Line 3: 



dw' (^^) 2?r(6 - 4/3 47 -f 3e - 2-f^)U[w'9, w']d,U[w'9, w'] 



Modulo the coefficients 7, /3, and e then this term scales like {rH/X)U'^, which is of order 10^^ for the dominant 
mode. An alternative theory of gravity would have to produce PPN-like parameters of magnitude 10"* or greater 
in order for this term to be as important as the Born correction. On the other hand, if the PPN-like parameters 
are < 10^ always, then this term may be difficult to measure except on very large scales [H] where w ~ A and 
the Born correction is suppressed by comparison. 



Line 4: 



j^'dw' (i^) + j)u[w'e,w'] j;;"' dw" d,u[w"e,w"] 



Noting that U ^ U/rn, this term scales like {rH/X)U'^, which is roughly of order 10 * for the dominant mode 
and even smaller for larger-wavelength modes. 



Line 5: J^' dw' (^) ^(1 + 27 + j^)U[w'e, w'] J^' dw" d^U[w"9, w" 



This term also scales like (th/X)!/'^. 
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Line 6: /J" dw' (^^) ^ {n x (v x V[w'e, w']^ } 



This is the gravitoniagnctic term. Note that W ■ V ^ If in the matter equations, meaning that \V\ ~ {X/rH)U. 
Therefore this term scales like U ~ 10~^, making it ~ 1% the size of the leading order term and 10 times smaller 
than the Born correction. 



Line 7: J^'' dw' 



W — w' \ 1 ( ^1 Ql 



^{e' - f3')d^^2[w'e,w'] 



Noting that $2 ~ t^^i this term scales like {rH/X)U'^ ~ 10 which is a factor of 10^ smaller than the leading 
order term. 

• Line 8: dw' (^^) ^,V,[w'lw'] 

Noting that \V\ ~ \V\/rH, this term scales like \V\ {\/rH)U. 

Therefore we can safely neglect terms 4, 5, and 7. We also neglect terms containing a time derivative of the PPN-like 
parameters, as arising from lines 4 and 8, which are suppressed by a factor {X/rn)- This leads to the final answer for 
the deflection angle, 



(GM)/^ 



where 



_ (I + 7) /• 



w — w 
w 



diU[w'0,w' 



is the piece that is first order in U. 

af\e,w) = / dw' 
Jo 
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(6-4/3 + 47 + 3e-272)^^^^ , ^ . 



2c4 



u[w'e,w']diU[w'e,w' 



(1 + 7)^ 



„4 ""-l 
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u[w'e,w'] 

Jo 



dw" {w' -w")djU[w"0,w"] 



is the piece that is second order in U, and 



af^\9,w)= r dw' f'^L^\ [Ijnx fv x •(/[«;' | - \Vi[w'0,w'] 

Jo \ W J \_C I L ] J i C 



is the piece that depends on the gravitomagnetic potential. 



D. The distortion tensor and convergence 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



The distortion tensor is 



where we note that 



Therefore the linear order piece is 



de. 



-U[w'e,w'] = w'djU[w'0,w'] . 



(1) ^ (I + 7) 

where we have defined the weighting function to be 



'i'ly = ^^V^ / dw'g{w,w')dijU[w'0,w'] , 
c Jo 



g{w,w') 



w'{w — w') 
w 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



7 



The second order piece of the distortion tensor is 

(2) _ (6-4/3 + 47 + 3£~27^) 



2c4 

(1 + 7)' 



dw' g{w, w') 



dw' g{w,w') diUlw'Ojw'jdjUlw'Ojw'] + U[w'9,w']dijU[w'9,w' 
„g{w',w") 



dw" 



w'd,jkU[w'9, w']dkU[w"9, w"] + w"d,kU[w'9, w']djkU[w"9, w' 

(2.24) 



and the gravitomagnetic piece is 



(GM) 



dw'g{w, w') 



\(rix [v X V\w'9,w']\ ] - M9,w'] 
r L L J J Li c** ' 



The lensing convergence of the hght rays from the redshift shce corresponding to the distance w is 



1 



(2)/ 



.(GM)/ 



containing the foUowing pieces: 



at first order in U , 



dw' g{w,w')\j'^U[w'9,w'] 



(2.25) 



(2.26) 



(2.27) 



K^^\9) 



(6-4/3 + 47 + 3e-272) 
(1+7)' /•» 



dw'g(w,w') / dw 



dw'g{w, w') I \VU[w'9, w']J + U[w'9, w']\/'^U[w'9, w'] 
g{w',w") 



w'dukU[w'9, w']dkU[w"9, w"] + w"d,kU[w'9, w']d,kU[w"9, w" 



(2.28) 



at second order in ?7, and 



dw'g{w, w')V 



-T,n X 



1 A 



VxV[w'9,w'] ~—V[w'9,w'] 



(2.29) 



to first order in the gravitomagnetic potential Eqs. ( 2.26|2.29 1 comprise the first main resuh of our investigation. 

All of these pieces of the convergence depend on cosmological analogues of the PPN parameters. Noting the 
dependences, we expect three-point statistics of the convergence field to generically be functions of the parameters 
appearing in k*^^' and k^^-* above. As an example, consider one of the simplest statistics, the skewness of the convergence 
field i22|, 5*3 = , where the angle brackets denote averages over the survey area. To lowest order, the 

denominator depends on the first order convergence given by Eq. ( 2.27[ ), and numerator depends on a combination 
of the first and second order pieces: (k^) = 3 (k^-'^^k^^^k^^^) + higher-order terms. Hence, the 2PN lensing deflection, 
including the cosmological analogues of the PPN parameters /3 and e, contributes to the three-point function of the 
convergence. 



III. THEORIES OF GRAVITY 



The indications of new gravitational phenomena are manifest in weak lensing through the cosmological analogues of 
the PPN parameters 7 at linear order, and (5 and e at second order. We have shown that lensing depends in particular 
on the combination (6 — 4/3 + 47 + Se — 27^). It turns out that, for scalar-tensor theories of gravity at least, the 
combination 3e — 4/3 depends simply upon 7. That e does not introduce a new degree of freedom has previously been 
noted [23]. We will show that the dependence of e upon /3 is precisely cancelled when the sum of metric coefficients 
5oo + 9xx required for lensing is taken. This means 7, or its cosmic analogue occasionally referred to as "gravitational 
slip," is all that is required to model weak gravitational lensing to second order. This applies to scalar-tensor theories, 
including specific theories such as the well-studied Jordan-Brans-Dicke (e.g. Ref. |TT]), /(i?) [inj, and chameleon or 
"symmetron" [24 theories. 
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Consider a scalar-tensor theory of gravity with the action 

1 ( m 



167r 



(t)R 



r 



(3.1) 



Jordan-Brans-Dicke theory corresponds to the case A = 0. The f{R) theory corresponds to setting i? = 0, A = 
^R{dhiR/d\Ti f — 1), and (j) = f^R. The symmetron model corresponds to a Higgs-like potential for V{(t)) where 
A(0) = —(j)V{(j)). All theories yield the same dependence of e upon (5 and 7. We now focus on the f{R) theory. 



A. f(R) Theory and Field Equations 



Consider the action 

for which the field equations are 

FRfiu ~ 
UF = 



d x\J —g 



IGtt 



(3.2) 
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1 /Stt 
3 I c 



F. 



{f-RF)g^, 



^T + 2f~RF 



(3.3) 



where F = df /dR. The second equation, equivalent to the scalar field equation of motion, guides the evolution of 
a new degree of freedom of the gravitational field as the Ricci scalar curvature "comes to life" as a dynamical field. 
Since the 2PN expansion in this theory has not been widely explored, we will provide some details of the calculations. 

We consider a spacetime filled with non-relativistic matter in order to model the case of a static, spherically- 
symmetric mass M. Hence, the stress-energy tensor is 



T, 



(3.4) 



The fluid four- velocity satisfies = — 1 in the full spacetime. We define the "conserved mass density" p* = py/~gu'^; 
see Ref. TT], Eq. (4.77). By ignoring the cosmic fluid, consistency of the unperturbed field equations requires 
f{R = 0) = 0. A full treatment including the cosmic fluid has been carried out elsewhere (see also [?7H25] ). Since 
our focus is on scales within the Hubble horizon, as argued in the previous section, then for clarity we can simplify our 
calculations and use a static spacetime metric. Note that we have also carried out an analysis in a Robertson- Walker 
background, which agrees with the simpler results presented below when we restrict our attention to length scales 
within the Hubble horizon. 



B. Metric to order 0(1 /c^) 
We now write the spacetime metric to order 0{l/c^) in the form 

ds^ = -(l- \Ai] c^df + ] dx^ 



(3.5) 



Hereafter the subscripts indicate the expansion order in inverse-powers of . The field equations give 



IGtt 



V^V-'(yli - 2Bi) 



/(O) 
3/"(0) 



3/"(0) 



P 



(3.6) 



where the Ricci scalar curvature to this order is i?i = V^(Ai — 2Bi). The exterior solutions for Ai and Bi are no 
longer simply proportional to 1/r outside the mass source; the loss of Birkhoff's theorem means that we must take care 
to match the exterior solution onto the interior solution. See Ref. [25 for a detailed exposition of this theory to order 
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0{\/c^). To be definite, let us consider a constant-density source of radius tq where p* = 3Af/(47rrQ). Furthermore, 
let us define = /'(0)/(3/"(0)). Now we can show that for r > tq 



4M 
/'(0)r 

4Af 
/'(0)ro 



^l(mro)fco('7^r') 



(3.7) 
(3.8) 



where 1^(2;) and fc£(x) are modified spherical Bessel functions of the first and second kind. Solving for the desired 
potentials, we find 



^1 = 77 



4M 



/'(0)r 
4Af 



1 H ii(mro)ko{mr) 



./'(0)r 



-ii(mrQ)kQ{mr) 



(3.9) 
(3.10) 



So if we equate Ai = 2U and i?i = -fAi then we infer 



G- 



1 



(l+7)/'(0) 



1 - 



2^ii{mro)ko{mr) 
- + 7-ii{mrQ)ko{mr) 



(3.11) 



In the limiting case mr ^ 1, ko^mr) — >■ and we recover GR with 7 = 1. In the case that mr, mro ^ 1, we find 
7 = 1/2. Moreover, 7 grows monotonically from 1/2 to 1 as r grows from r <^ to r ^ m~^. This highlights 
a feature of /(i?) gravity, that 1/2 < 7 < 1. We also learn that the contribution to the curvature does not vanish 
outside the massive body: it decays like a Yukawa potential with mass m. This spatial dependence of 7 manifests 
itself as small-scale fiuctuations which are below the resolution of the particular weak lensing observations of interest 
here. 

We also notice that the sum of metric pote ntial s that appear in the lensing equations such as Eq. (2.8 1, i.e. Ai+Bi = 
2(1 -I- 7) [/ to order ©(l/c^) as given in Eq. (3.7 1, decays like as in GR. Over the range of length scales for which 



these equations are valid, this means gravitational lensing is insensitive to the Yukawa potential - lensing does not 
feel a fifth force. On the other hand, estimates of the mass AI based on orbital dynamics rely on Ai and do feel the 



Yukawa potential in Eq. (3.9). 



C. Metric to order 0(l/c^) 



To order 0{l/c^), the spacetime metric gains a term: 



1 



Ai c^dt 



1 



Bi dx^ 



The off-diagonal, t — i component of the field equations yields 



167r 



P Vi 



rVicdt dx^ 



(3.12) 



(3.13) 



where Vi is the leading contribution to the spatial component of the source 4- velocity, Ui — Vi/c + O^l/c"). The 
solution is Vi — — 4Vi//'(0), where we define Vi implicitly as the solution to V^Vi — —A'Kp*Vi. Substituting for /'(O), 
we find the familiar result that Vi = —2(1 -f ^)GVi. The appearance of 7 is somewhat deceiving, however, since it is 



cancelled out by G, cf. Eq. (3.11). Consequently, the vector potential is also insensitive to the Yukawa potential in 
Ai. 



D. Metric to order 0(1/0*) 

As we proceed to higher order, we will ignore the time-dependence of the metric potentials and also ignore the 
contribution of the source velocities. What this means is that our results up to order 0{l/c^) are generally valid for 
a non-relativistic source, but higher order results are only valid for a static source. 
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We write the spacetime metric to order 0{l/c^) as 



1 - -^Ai - - 



A2 I c^dt'^ + —Vicdtdx' + ( 1 



(3.14) 



For our purposes it is sufficient to examine only the time-time (/i = = 0) component of the perturbed field equations, 
whereupon 



V2(A2 + B2) = (aI + A,B, - ^Bf + -\a,R, 

4 \ 2 



1 

61m? 



Ri 



2n 



P* [Ai+ 2Bi 



3m- 



:Rl 



(3.15) 



Defining $2 implicitly as the solution to V^<i>2 = — 47rp*J7, then the solution to Eq. (3.15) is 

1 



^2+^2 

(/?' + £')^2 



, ,Al+ AiBi - ^Bf + -\AiR, )+(/?' + s')^2 
4 V 2 m-^ 



1 



/'(O) 



/'(O) ^2 



247rTO2 ^ 



(3.16) 



Again, this is the sum of metric potentials we require for gravitational lensing. Since Ri/m? = |(1 ~ 7)^1 1 '^e find 



A2 + B2^ — ij(5 — 7 — 37^)?/^ + {P' + e')<^2- Comparing with the general metric Eq. (2.1), we find 

3£ - 4/3 = -5 + 7 + 37^ . 



(3.17) 



This combination, 3e — 4/3, is what appears in the lensing convergence, Eq. (2.281. Hence, in the context of this theory 
of gravity, the 2PN lensing terms depend solely on 7. 

Phenomena other than lensing depend upon /3 and e separately, so we turn to the spatial components of the 
gravitational field equations to find 



V^B2 



1 



V^V^(A2 - 2^2) . ^^p^B. + l^^Bl) - ^V^(i.?) + ^C. 



r(o) , 

/'(O) 



1 



(3.18) 



where 



C2 = -V2v2(yl? + AiBi + 355^) + ^/'^{AiV^Bi - 2Ai\l'^Ai + 4SiV^Ai - SBiV^Si) 
-^(V^AiV^Bi - V^AiV^Ai + 2V^BiV'^Bi) . 



(3.19) 



In the GR limit, takin g m — » cx) and /'"(O) — ^ 0, the first term on the left-hand side and the first two terms on the 



right-hand side of Eq. (3.18) dominate. As a result, B2 = iU — $2 and A2 = —2U — 2$2, meaning that /3 = e = 1, 



as expected. In the limit m 0, and for r^/"(0) ^ /'"(O), then the 2nd term on the left-hand side and the fourth 
term on the right-hand side of Eq. (3.18 1 dominate, leading to /3 — 1 and e ~ 1/12. The combination 3e — 4/3 declines 



monotonically from —1 for GR down to —15/4 for 7 = 1/2 



E. Other theories of gravity 

The f{R) theory predicts 1/2 < 7 < 1. However, let us briefiy consider a gravitational theory described by the 
action Eq. ( |3.1[ ) with a potential A(0) and kinetic coupling 6'((/)). We presume that A, 0, and </> can be expanded 
perturbatively around their values on the backgro und (j ) — 0Oj as described in explicit detail in Ref. [26) . Next, 

the equations to order 0{l/c^) are 



HIE 



following the same procedure as carried out in Sec. 

V2(Al+Bi) 

V^Ai - Bi) - m\Ai - B,) 



167r 

<l>o ' 
167r 



(bo 3 + 20(00 ) 



(3.20) 



where = —2(f>QX"{(j)o)/[3 + 29{(j)o)]. Considering a constant-density sphere of radius rg in the limit rg 0, the 
solution for 7 is 



7-1 



l + e""-[3 + 26l(0o)] 



(3.21) 
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We assume that (j>o\((j)o), 0qA'((/)o), and the cosmic energy density are all negligible in comparison to the matter 
source p. Consequently, the only way to achieve a large value of 7 somewhere while still satisfying the Solar System 
constraint 7 = 1 + (2.1 ± 2.3) x 10"^ [501 131] is if -2 < ^'(^o) < -3/2 and mr > 1 on the scale of the Solar System. 
But even in such a case, at larger radii beyond the Solar System, though smaller than the Hubble scale, 7 asymptotes 
even closer to the GR value. 



We can again see, in Eq. (3.20), that the lensing potential is insensitive to the Yukawa potential or fifth force, and 
instead decays like 1/r as in GR. What lensing does feel is the conformal factor (p^ that modulates the strength of 
the source in the above Poisson-like equation. Again, if that conformal factor can be determined elsewhere, such as 
by measurements of dynamics at radii within the fifth-force cutoff range, then lensing can give an undistorted view 
of the mass despite possible hindrance by a possible fifth-force. 

To extend our study, we may consider the chameleon mechanism, whereby the non-linear solution to the OF equation 



(3.3) yields agreement with GR in high density regions such as the Solar System and Milky Way, but departure from 
GR leading to cosmic acceleration in lower density environments on larger scales. The requisite equations are little 
different from what we have already derived, but necessitate that we keep certain background terms that we have 
ignored in the preceding analysis. However we can borrow the results of Ref. f35^ (hereafter HS) to illustrate our 
points. To translate between the different notations, Ai = 2Ahs — "^IShs, Bi = 2AhS: and f = {R + Jhs)- For a 
Universe containing only non-relativistic matter, the sum of potentials that appear in the lensing equation satisfy the 
equation 



2 , . ^ . IStt 1 



V^(Ai+i?i) = --^p+-A , (3.22) 



whereas the individual potentials satisfy 



where A = R — f{R)/f'{R). Here R is the solution obtained from solving Eq. (3.3) in the environment of stars, 
galaxies, and clusters as described by the non-relativistic mass density p. In order to satisfy Solar System constraints, 
/ ~ R/G is required to lock the curvature to the density. To borrow the language of HS, this criterion helps define 
a critical value of curvature, above which GR is recovered. In such regions, A <^ R, Sirp/ f (R). In viable models of 
f{R) gravity, then the proportionality between the lensing potentials and the mass distribution is simple l//'(i?). In 
regions where 7 = 1/2, leading to Ai — 2i?i, then we can see that A, i? <C 8TTp/f'{R). This means the sum of the 
potentials is still generated by the density p and not the novel curvature via R, A. Furthermore, if the bulk of the 
mass is located in a region in which 7 = 1, then the gravitational constant serving as the constant of proportionality 
between the lensing potential and mass is the same as that safely measured on smaller scales. 



The relationship between (3 and e given in Eq. (3.17) is valid for scalar-tensor theories with a single scalar field, 
as has been shown in Ref. [23^. For a single scalar, there is a single new free function whose effect on the metric 
potentials is completely described by a single parametric function at the order C(l/c^), namely 7. In the presence 
of additional scalars that couple non- minimally to gravity (e.g. |33j). new functions may be necessary. It appears to 
be an open question whether vector-tensor and massive gravity theories, such as the galileon |34| . will also require 
additional functions to describe the departure from GR. Interesting signatures of galileon theories have recently been 
found for cluster lensing [55] . 

IV. DISCUSSION 

We have explored what happens to the weak lensing distortion field at second order when the nonlinear behavior of 
gravity is altered. If we parameterize such deviations in a cosmological adaptation of the PPN framework, i.e. starting 



from the metric (2.1), then we find corrections to the nonlinear weak lensing convergence that depend on the PPN 
nonlinearity parameters /3 and s. However, for a wide class of scalar-tensor gravity theories, we find that this 
combination generally cancels, thereby leaving the nonlinear correction to only depend on the "gravitational slip" 7. 
We have argued that the time-evolution of the PPN-like parameters leads to a negligible correction on scales well 
inside the Hubble horizon. We further derive the spatial dependence of 7, which is expected to be important only on 
sub-galactic scales and therefore would not impact the weak lensing observables di scus sed here. 



Our results can be cast in terms of an effective A/nl- Returning to the metric (1.31, then the sum of the 0{l/c^) 



potentials as appears in the lensing deflection is —2f3U + ^eU = —2(1 — A/nlJC/^ + 2(4 -I- A/nl)C^ • Matching 
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this sum with A2 + B2 as computed in Section HID, we determine that the additional non-Gaussianity introduced by 
gravity beyond GR is 

A/nl = ^(1 - 4^ + 3e) = ^(37' + 7 - 4) . (4.1) 

For the range | < 7 < 1, we find |A/nl| < 1- The tightest constraints to date come from the CMB [S^: /nl — 32±21. 
Based on the scalings of the various contributions to the deflection angle given in Sec. |II C[ a scalar-tensor theory 
would have to predict 7 ^ 10^ for these effects to be observable. But such a large value of 7 is not possible, and indeed 
this parameter has been constrained to be close to the GR value (e.g. [9]). If this theory satisfies Solar System tests 
of GR then the non-Gaussianity of gravitational lensing relative to GR on scales below the Hubble scale is too weak 
to observe. However, it is worth noting that modifications to nonlinear gravitational dynamics can further impact 
/nl via modifications to nonlinear structure formation. Overall, our results suggest that looking for the signs of new 
gravitational physics in the linear regime, in the evolution of large-scale clustering over cosmic time scales, e.g. as a 
consequence of "gravitational slip" , remains a valid strategy. 

If the Solar System tests are satisfied, then the gravitational constant measured in the Solar System is G = 1//'. 
This is the same gravitational coupling as appears in gravitational lensing: for a mass M, the deflection relies on the 
sum of potentials Ai -\- Bi ~ 4Af//'r. This means that sub-horizon gravitational lensing may provide a clean, direct 
method to measure mass in spite of a potential fifth-force resulting from a scalar-tensor modification of gravity. 
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